setwd("~/Cline_disp_evol/Balkau_orig/ran_outs")

concat<-c()
for (ii in 1:107){
  load(paste("balkau_orig", ii, ".RData", sep=""))
  concat<-rbind(concat, ptable)
  rm(ptable)
}

concat<-as.data.frame(concat)
concat$mdiff<-concat$m-concat$m.star

png(filename="../Figures/diff.png", width=12, height=12, units="cm", res=300)
plot(concat$mdiff, concat[[5]], col=terrain.colors(n=4)[factor(concat$m)])
dev.off()